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ABSTRACT 


Recent research which contributes to the prediction of the acoustic impedance 
of duct liners is reviewed. This review includes the linear and nonlinear 
properties of sheet and bulk type materials and methods for the measurement 
of these properties. It also includes the effect of gracing flow on the 
acoustic properties of materials. Methods for predicting the properties of 
single or multilayered, point reacting or extended reaction, and flat or 
curved liners are discussed. Based on this review, methods for predicting 
the properties of the duct liners which are typically used in aircraft engines 
are recommended. Some areas of needed research are discussed briefly. 
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SUMMARY 


The methods of predicting the acoustic properties of duct liners is reviewed. 

The review is divided into a discussion of typical acoustic materials aftd a 
discussion of methods of predicting the liner properties from these material 
properties. Acoustic materials are usually described by their acoustic 
impedance. For thin sheets of materials, this impedance relates the pres- 
sure differential across the material to the velocity through the material. 

For bulk materials, the impedance relates the pressure gradient to the 
velocity. There are two commonly used methods of measuring impedance. 
These are the standing wave method and the two microphone method. Acoustic 
materials exhibit significant amounts of nonlinear behavior at the high inten- 
sities usually found in aircraft engine ducts: however, this nonlinearity iS 
well described by the measured flow resistance of the material. Also, 
tangential flow has a large effect on material properties. The best available 
information indicates that these effects can be accounted for by the material 
flow resistance in the presence of tangentialflow. In most practical problems, 
it is necessary to linerize the material properties by use of an rms velocity 
in the material, however, this velocity must be found by an iteration procedure. 
Bulk materials can be characterized by two measured parameters, the imped- 
ance and wave number. The properties of complex duct liner configurations, 
are predictable once the liner material acoustic properties are known. 
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SYMBOLS 

acoustic resistance constant 
‘‘/Oa'a 

nonlinear acoustic resistance constant; liner surface radius 
in Equation 46 

b*/p ; kb* 

cl a 

orifice discharge coefficient 

speed of sound in ambient or duct medium 

distance from face of material to first minimum of Lp in standing 
wave measurement method 

distance between minima of Lp in standing wave measurement 
method 


liner backing cavity depth or layer thickness 


hole diameter of perforated sheet 


k d 


a 


* 

h 


frequency 


material acoustic reactance constant 




woi' 


a 

h* 

b 


^ i 


nonlin6ar itiatcrial acoustic reactance constant 


h*/p. 


K 


"U) 


CO 


rm 


zm 


Lp 


M 


Mw 




grazing flow constant 


wave number 


k*/k. 


= u)/c. 


complex wave number constant for bulk materials 


ki/k^ 


wave number in extended reaction liners 


acoustic mode axial wave number in duct 


standing wave ratio, dB 


sound pressure level, dB re 2 x 10 N/m 


f'"' thickness of porous or perforated sheet material 


k f 
a 


Mach number of mean tangential, grazing or duct flow 


material acoustic inertance 
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m 


circumferential mode number 


N 


number of sheet material layers in a multi-layer liner 


number of frequency bands 


p’-‘(t) 


p(t) 


acoustic pressure 


P 

P 


* 


mean static pressure 


*/ 2 


'U> 


spectral (complex) acoustic pressure 


# / 2 
Pc. /Pa<=a 


nth harmonic of , where cj = nt.O| 


- OU) 


incident spectral (complex) acoustic pressure 


• OW 


P /pc 
^ ow'^a a 




grazing flow pressure fluctuations 


R 


steady flow resistance normalized by p c 

a a 


temporal (instantaneous) acoustic resistance 


R. 




R 


va) 


sheet material acoustic resistance = Re(AZ^ ) 


R 




IX 


R 


R 


Bheet material specific acoustic resistance at nth harmonic r.« 

sheet material specific acoustic resistance basecj on steady flow 
resistance, R 


«h 

r 

e 

t 

t 

u 

v(t) 


V 


specific orifice or hole acoustic resistance 

effective honeycomb cell radius 

time 


= Wt 


mean tangential, grazing or duct flow velocity 
acoustic particle velocity normalized by c 

a 

steady-flow Mach number normal to liner 
spectral (complex) acoustic particle velocity 


u> 


/c 


a. 


n 


>2' 

xt 


nth harmonic of where w = nwj^ 

total effective particle velocity 

temporal (instantaneous) acoustic reactance 


X-; 


sheet material acoustic reactance = -Im(AZy ) 


“w 


x„/P^S 


X 


n 


X 


sheet material specific acoustic reactance at nth harmonic nw 

measured and/or calculated sheet material specific acoustic 
reactance 


1 


X 


axial coordinate in standing wave tube 




0 ^ 


X 


AZ 


^ Z 
iJ 

z; 
z< 

tk z 

n 


p 

t 

' (a) 




sheet material acoustic impedance = - IX^ 

Z^Z*/p^c^ zK^-iX^ 

instantaneous acoustic impedance of sheet materials 

liner acoustic impedance = AZ* + zt,^ for a single layer liner 

</PaS 

sheet material specific acoustic impedance at nth harmonic nuj 


^ U) , 


acoustic impedance of single- layer cavity 




/ 

^«c^^Va 




spectral (complex) distributed impedance of a bulk acoustic 
material 


w 




attenuation factor 


S|* 


p 

p 

6* 

6 


acoustic admittance 


❖ 


p c P 

^a a 


linear boundary layer thickness 


k 6 
a 




boundary layer displacement thickness 


k 5* 

'^ar I 


xi 


X 


6ound wavelength 


(J 




U 


kinematic viecoeity 
density of ambient medium 
sheet material effective density 
= Pe/Pa 

open area of perforated sheets porosity of material 
circular frequency 
fundamental frequency 
reference frequency 


Sub sc ripts 

i denotes frequency band i (in place of w' 

4 - refers to steady flow suction normal to the liner 

(i. e. , from duct into liner) 

- refers to steady flow blowing normal to the liner 

(i. e. , from liner into duct) 

Supersc ripts 

refers to nth lay» r in a multilayer liner 

* denotes a dimensional quantity: exceptions arc the common 

symbols c^, k^, 


xii 


Operators 


Re, Im real, imaginary part of 

^ difference in a quantity 

^ gradient vector 

2 

V divergence (gradient) or Laplacian 

3/ 9t partial differentiation with respect to time, t 

^ temporal (instantaneous) impedance operator 

R* + X* 8/8t* 

^ ^t 
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Section 1 
INTRODUCTION 

Duct liners are used to reduce noise in a number of situations. One 
application for duct liners is the nacelles of modern aircraft engines. In 
this environment, the liners are exposed to high-intensity sound, (Refs. 3, 
2, ?), high-speed flow, and large temperature variations, B cause f non- 
acoustical requirements, only a few materials have ever been quali :i(;c‘ for 
service in aircraft engines. A perforated sheet is one material 'v.<’ch has 
been used, porous metals and rigid fiberglass sheets ■■ ■ . oth't 

One goal of aircraft noise research is to dev lop the ability to predict the 
acoustical performance of new nacelle duct treatments. This prediction 
ability must be based on the properties of duct liners; hence, the present 
paper considers the scbproblem of predicting the properties of duct liners. 
In any prediction scheme, experimental data must be included at some point, 
for example, the stiffness of a structure is predicted from the experimental 
knowledge of the modulus of elasticity from which the structure is made. 

Thus, part of any prediction scheme is the identification of those properties 
which are to oe measured. 

This report reviews the state of knowledge of the properties of acoustic 
materials and the way these properties are measured. The physical 
variables which affect the acoustical properties of materials are discussed 
and measurement methods are recommended which account for the presence 
of the most important physical effects. Then, methods for predicting the 
acoustical properties of complex duct liner configurations are given. 
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Section 2 

ACOUSTIC MATERIALS 

The propagation and attenuation of sound in lifted flow ducts is understood in 
physical terms such that a mathematical description falls naturally into two 
parts. One is generally concerned with how sound that is generated within 
the duct propagates through the duct flo\tr field and out through the duct term- 
ination into the far field, while the other is centered on the absorption of 
sound by the duct boundaries. Both areas have been studied intensively 
although the former has attracted more attention recently since it has been 
realized that the sound radiation can be sensitive to mean flow conditions 
within and exterior to the duct even when dissipation processes within the 
flow are ignored. The absorption or dissipation of acoustic energy by the 
lined duct surfaces is historically a well established subject, but it has been 
known for some time that in the present context the classical models of 
sound absorption must be considerably modified to take i.ito account high- 
intensity sound and grazing flow effects. 

In this section ',;e are concerned with the way in which the physical process 
of sound absorption by an acoustic material is described through the imped- 
ance concept. In general this concept allows us to link the absorption and 
propagation processes without the need to consider how sound absorption is 
related to the physical properties of the liner materials (e. g. , porosity, 
surface density, fiber size). That is then a separate problem which is 
outside the scope of the present effort, but nevertheless is one of consider- 
able importance. We also review the largely empirical methods which 
attempt to represent the observed influence of sound intensity and grazing 
flow on the sound absorption processes through a modification of the imped- 
ance parameter. In a sense these real effects affect the convenient division 
of the mathematical description into two parts; that is, in practice, the 
absorption process cannot be specified until the local sound field is known 
but this depends on the propagation processes (as well as the local sound 
field) so that one process depends upon the other. 
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2. 1 LINEAR IMPEDANCE 

Impedance is generally defined in acoustics as the latio of the acoustic 
pressure and the acoustic particle velocity at a particular frequency, (in 
the following it is implicit that the impedance is always a function of fre- 
quency), hence, it is in effect a vector quantity and a general property of 
the sound field anywhere in the duct or in the radiation field. As in linear 
electrical circuit theory it is a useful parameter in acoustics since the 
dynamics of these systems is, by definition, independent of absolute ampli- 
tudes and only the ratios of the variables are important. 


2. 1. 1 Sheet Materials 

The acoustic impedance of a material is defined in Reference 4 as the 
complex ratio of sound pressure to the component of particle velocity normal 
to the surface of the material. 


Ui 





( 1 ) 


In equation 1, both p* and are understood to be the amplitude of complex 
harmonic motions in time with the time factor exp (- iwt). Reference 4 
describes the classical impedance tube method of measuring the impedance 
of acoustic materials. 

Equation 1 gives the total impedance at the surface of a material, such as 
in Figure la. 

For thin sheets of material, it is convenient to consider the sheet to have 
zero thickness and define the impedance by 

A 7* - ^ ~ ^ Pm (2) 

^ - if ~ 
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W -:*■ MATERIAL 


A. SULK MATERIAL SAMPLE 


P„M 


CAVITY 


B. SHEET MATERIAL SAMPLE 

Fisum 1. ImiMdinn Conctpti for Bulk and ShMt Matarialt 

This definition amounts to subtracting out the impedance of the backing cavity 
in Figure lb. This backing impedance is approximately i cot (k* d*), where d* 
is the backing depth. 

At sufficiently low sound pressure levels the impedance of a thin sheet is of 
the form 




where .M the material inertance, is usually a fairly weak function of 

'i' 

\ * both and X^^are a function of the niaterial geometry and the 
mean properties of the local fluid, e. g. , its coefficient of viscosity and mean 
density as well as frequency. Thus, although can be measured quite 

easily, data should be corrected if fluid properties in the measurement and 
application environments are significantly different. It is beyond the scope 
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of the present report to review the functional dependence of on the 
material geometry and the fluid's mean properties although it is given when 
we consider the special case of perforate materials below since it is of a 
particularly simple form. 

2.1.2 Bulk Materials 

Scott (Reference 5) has shown that the propagation of sound waves in iso 
tropic porous materials can be characterized by two complex constants. 

Scott used a complex wave number and density; however, other pairs of 
constants may be used also. For example, Bies (Reference 6, ch. 10) uses 
a complex propagation constant and characteristic impedance to characterize 
bulk materials. 

In this report, we will use a distributed impedance, or impedance per unit 
length, to relate the pressure gradient to the velocity in an isotropic porous 
medium. 

^ « * * 

= 0 (3) 

The form of equation 3 is used here because of the analogy to equation 2. 

The gradient of pressure is approximated by - A p^/f where t is the 
thickness of a thin sample such as in Figure lb. Equation 3, together with 
the wave equation given by Scott (Reference 5). 

_2 * . * 

^ P« + ko, = 0 (4) 

completely describe the properties of waves in infinite media, porous or 
otherwise. In the limit of 100 percent porosity the porous medium is the 
ambient medium and then 
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which is a real number^ and 


* 

= -iwp- 


which is a negative imaginary number. 


2.1.3 Measurement Methods 


The simplest method of measuring the impedance of a material is the 
standing wave method described in Reference 4. Figure 2 (from Reference 4) 
shows the apparatus for this measurement. 

This method is attractive because it does not require a lot of sophisticated 
equipment for its use. A speaker, driven by an audio oscillator (and perhaps 
an amplifier) sets up a standing wave pattern in the tube in which the sample 
is placed. The experimenter moves the microphone away from the sample 
face, noting the minimum sound pressure level Lp (Dj), the distance to that 
minimum Dj, the maximum level Lp (Dj) +L, and the distance between minima 
D2- From these measurements, and the assumption that the waves in the tube 
are of the form exp j±ik*x*} the impedance of the sample may be deduced. 
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One difficulty which occurs with the standing wave method is that the 
effective wave number k* is not real due to absorption of sound by the tube 
v ails. This situation was analyzed by Lippert (Reference 7) who considered 
waves in the tube of the form exp {+(lk*x* + a*x*)) where is an attenuation 
factor due to the tube wells. In this case, it is necessary to determine the 
envelope of the pressure minima* as shown in Figure 3, to correctly find 
the impedance of the material. 

Thin sheet materials are characterized by the impedance which is given by 
Equation 2. This impedance is found, using the standing wave method, by 
subtracting the impedance of the backing cavity (Figure lb) from the imped- 
ance of measured impedance. Since the impedance of the backing cavity is 
approximately i cot d, where d is the cavity depth, the measured impedance is 
large when d is n§ar mr, n = 1, 2. . . The measured impedance in front of the 
sample is consequently large and the sample inpedance, which is found by taking 
the difference between two large numbers, will not be accurate. The obvious 
solution to this problem is to hold d at ir/2 (the quarter-wave cavity) such that 
the cavity impedance is zero and the measured impedance is the material 

CMS 
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Impedance, While this is simple as a mathematical concept, it introduces 
the practical difficulty of adjusting the cavity depth to \/4 for each test 
frequency. 

A second method for measuring the impedance of a thin sheet is the 
two. microphone method. As the name implies, an additional microphone is 
added to th, test apparatus in Figure 2. This microphone is fixed at the end 
of the cavity and provides a measure of both the velocity through ;he thin 
sample and the pressure immediately behind the Sample. The method also 
requires an additional piece of equipment, a phase meter, to determine 

the complex ratio of the pressure in front of the sample to the pressure at 
the wall of the backing cavity. 

M.lling (Reference 8) ha. given an extensive discussion of the design and 

operation of a precision impedance tube using the standing wave method of 
measurement. 


MelUng also provides a brief discussion of the two- microphone method and 
note, that its principal difficulties are errors in phase measurement, which 
IS an instrumentation problem rather than a fundamental defect of the 
method, and errors due to placing one microphone at the face (in front of) 
the sample. Melling was studying the impedance of perforated sheets and 
the pressure field in the neighborhood of an array of orifices is very 
nonumform. This problem can be overcome by moving the microphone 
back from the sample face a fraction of a wavelength. When this is done 
the formula for the impedance of the sheet is slightly more complex than if 
the microphone were right at the sheet face, but the sheet impedance is 
>till given as a funcUon of the complex pressure ratio between the micro- 
P ones. Like the standing wave method, the two microphone method works 
best when the backing cavity is one-guarter wavelength in depth. 

Bulk materials must be characterised by two complex constants, as shown 
by Scott in Reference 5. There is a wide range of choices available for 
which two constants are used, such as compressibility, sound speed, wave 
number, characteristic impedance, or distributed impedance, however, 
t ere are only two independent constants. We have arbitrarily selected the 
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wave number (Equation 4) and the distributed impedance (Equation i) for use 
in this report. Since there are two independent complect constants which 
characterize the material, at least four physical quantities related lo the 
material sample must be measured. Scott (Reference 5) obtained the complex 
wave number (propagation constant) by moving a microphone through the 
bulk material while observing the signal of the microphone output and the 
signal from the driving oscillator on a dual beam oscilloscope. Figure 4 
shows a schematic of Scott's apparatus. The wave length was found by 
moving the microphone until 360® phase shift occurred, and the attenuation 
of the wave could be observed by the change of the microphone output as a 
function of distance-. These two constants are easily related to the complex 
Wave number of the material. A modern version of Scott's experiment 
would use a two microphone technique with a phase meter used to determine 
the wavelength. The second complex constant used by Scott to characterize 
the material was the characteristic impedance. This is the impedance at the 
front of an infinitely deep sample of the material which can be measured by 
the standing wave method. 








2. 2 NONLINEAR IMPEDANCE OF THIN SHEETS 

Up to this point it has been implied that the material impedance is independent 
of the absolute amplitudes Of the acoustic motion; by definition this should 
be so since the propagation processes are described by acoustic, or linear 
equations. However, it has been pointed out that the propagation processes are 
not considered within the material itself, as such; this is avoided by the use 
of the material impedance concept. It turns out that the fluid motion near 
and within many lining materials of practical interest, at the sound pressure 
amplitudes encountered in aero>engine ducts, is not a linear process and 
hence the impedance concept has to be used with care. 

An effective liner impedance needs to be redefined; this is outlined in this 
section and the necessary types of modification of the linear material 
impedance parameters are described. It also appears that, although less 
well documented and understood, the effective liner impedance is also a 
function of the local grazing mean flow velocity. 


When the fluid motion within a lining material is accurately described by 
linear equations it follows that a single acoustic wave incident on the liner 
surface will be partly absorbed and partly reflected, the frequency of the 
reflected wave being the same aS that of the incident wave. In addition the 
impedance seen by the local sound field is independent of the incident wave 
amplitude. In practice the impedance of most lining materials is dependent 
on the incident wave amplitude and furthermore the frequencies of the 
motion are not confined to the frequency of the incident wave. The existence 
of nonlinear effects is corroborated by the zero frequency or steady flow 
case where it is Well known that the mean pressure differential across a 
porous material at a sufficiently high Reynolds number is proportional to 
the square of the steady flow velocity. Based on this zero frequency case 
it has been postulated by Zorumski and Parrott (Reference 9) that the 
instantaneous pressure differential Ap(t) across a thin porous material is 
related to the total instantaneous fluid velocity v(t) normal to the sheet by 

Ap(t) = R^[v(t)] v(t) + [v(t)] 1^ ( 7 ) 
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where [v] is the measured steady flow resistance at the steady flow 
velocity V, and subscript t emphasises the dependence on time, t. Since 
at high Reynolds numbers 

( 8 ) 

where b is a constant; it follows that (neglecting [v(t)] for the present) 
Ap(t) = b v(t) jv(t) I ( 9 ) 


Suppose now that only one frequence, w^, of fluid motion is present, 

v(t) = Vj sin(t) ( 10 ) 

(here t = u>jt ) then 

^p(t) = bvj sin(t) I sin(t) j (H) 


or 


+ bvj sin^(t), o < t ■< IT 


Ap(t) = 


- bVj sin^(t), n < t S 2it 

Taking the Fourier transform of each side We get, following Zinn 
(Reference 10) and Melling (Reference ll) 


( 12 ) 


8b ,2 

^Pl = 3?^l <13) 

where Apj is the Fourier component of the fluctuating pressure differential 
at the fundamental frequency Wj. It follows that the impedance expressed as 
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a ratio of the pressure differential and the particle velocity at the funda- 
mental frequency is 


AP, 8b V. 

P - * * 

Vj " 3ir 


(14) 


At sufficiently low Reynolds numbers R [v] is independent of v 
R^v] = a 


(15) 


and then R^ = a 


( 16 ) 


In general, both linear and non-lineaz effects are present so that 


R 


1 


21 a -f 



(17) 


This oversimplified example serves to illustrate two points: (1) that given a 
particle velocity motion at particular frequency, fluctuating pressures at 
other frequencies are generated when nonlinear effects are present (given 
by the Fourier components at higher harmonics of the expression forAp(t)) 
and (2) that an impedance can still be defined in the nonlinear case but it 
is proportional to the particle velocity amplitude. In general, the particle 
velocity motion will contain a wide range of frequencies and hence the 
impedance at one frequency is influenced by motion at all other frequencies. 
Thus, although the impedance parameter c^n still be defined in the same 
way, it will no longer be a function of the material properties alone but is 
a nonlinear function of the spectral distribution and level of the local sound 
field (References 9, 12). One other point that is obscured in these intro- 
ductory notes is that the particle velocity is hot known, a priori, rather it 
is an estimate of the sound pressure level. Thus, an iterative procedure is 
required in general, although with certain assumptions this can be replaced 
by the solution to a simple quadratic equation. 
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2.2.1 Peffoyated Sheets 

A common approximation for the Impedance of perforates Is 
= [a + bvj - l[g - hv^] 

where v , Is an estimate of the total effective particle velocity normal to 
the material surface; b and h depend on the material geometry and mean 
properties of the local fiuld but h Is also a function of v^. For all tain 
porous materials the parameter b can be determined from the steady flow 
resistance at large Velocity v. For perforate materials three similar 
empirical expressions for b are available, one is given by Mariano 
(Reference 13) as 


b = 1. 14 


, 2 -‘/ 


exp 


..5071 /dj^ 


1. 7 


•X- 

( 4 ) 


<r 

V 


(19) 


The other expression, given by Rice (Reference 14), is bound up with the 
definition of the total effective particle velocity which he gives as 



l i=l 


Vij= |vj| exp 

and V. is the effective particle velocity at frequency band i. With those 
definitions Rice (Reference 14) defines the parameter b as 

w I (22) 

There is a similarity between the two expressions if we ignore the exponential 
terms and note that usually « 1; then they differ by a factor 1. l4/a°- . 

Also the exponential term in Rice’s (Reference 14) expression has an argu- 
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ment which is of the same form as that in Mariano's (Reference 13) that is, 
the ratio of a material length scale (hole diameter or sheet thickness) and * 
the particle displacement magnitude ( |v*|/q). Mariano's (Reference 13) 
expression also contains another exponential term, which is a function of 
the hole aspect ratio f»!«/d*, which is absent from Rico's (Reference 14) 
expression. The explanation for the former term is given by Rice 
(Reference 14); he argues that when the fluid displacement is small compared 
with the hole diameter the "jet will not fully form" and hence the dissipation 
will decrease. The choice of the hole depth as the significant length scale 
is seconded by Melling (Reference 11) who has shown that the nonlinear 
resistance depends on the orifice discharge coefficient which, in turn, 
depends on the ratio of thickness to hole diameter. The origin of the terms 


in these expressions is, of course, in the well established, steady flow 
expressions. 


A link between these two expressions is found in that given by Armstrong 
(Reference 15): 


" 2</2 ‘-O^Sl^^jexpj^-.SOTZ (23) 

that is, it is very similar to Mariano's expression except for the exponential 
term representing the displaceme..* effect which is identical in form to that 
given by Wee (Reference 14). To summarise, there is clearly agreement 
that the parameter b should be proportional to the inverse square of the open 
area ratio but that is all, apart from the consensus that a correction should 
be appUed in the form ♦(ff./v*) where (* is a representative length or 
effective thickness of the material (Reference 16). More work is needeo 
to establish a generally accepted, precise form for this function and any 
Other corrections which appear to be necessary. 

The correction for the nonlinear reactance h is based on experimental 
observation such as by Ingard and Ising (Reference 18^ Zorumski and 
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Parrott (Reference % and MelUng (Reference 11), but no formulae appear 
to be generally available for this effec*'. 

The calculation of the effective particle velocity ifl, according t»» fioiv.e 
authors, bound up with the gracing flow effects. 

2.2.2 Uniform Poroua Material s 

At high sound pressure levels when nonlinear effects occur Zorumoki and 
Parrott (Reference 9) have verified experimentally that the basic assumpt;.'; 
of Equation 7 is correct, that is, the material properties are simply 
described in the time domain by the measured functions R^Cv] and [v]. 
These functions are the generalizations of the resistance and reactance user, 
to characterize the material in the frequency domain. The spectral 
impedance is not uniquely defined since it will depend on the incident 
pressure spectrum in a complicated way. For thin porous sheets of 
material, it is believed adequate to assume constant 

reactance measured in a low intensity test. The resistance [v(t)] can be 
measured in a flow resistance test since is a function of velocity only. 

Assuming a line spectrum where Ap(t) and v(t) are given by 


V 

n 

an approximation to the spectral impedance is 

where v (t) is an estimated velocity. The velocity estimate may be the 
result of solving the problem with the assumption that is a constant which 

is, in turn, estimated. 

Very little is known about the dependence of the resistance (which is 
also the flow resistance) on more fundamental properties of porous 
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materials. Zwicker and Kosten (Reference 17) and Bies (Referetice 6) 
express the flow resistance in terms of such things as porosity, bulk 
density, and structure factor, but these variables must be measured so this 
equation is just pushing the measurement problem from one spot to another. 
It seems preferable to regard R^ and as fundamental properties of and 
acoustic material and fa:e the problem of measuring these functions 
directly. 

2.2.3 Measurement Methods 

Measurement methods for nonlinear thin sheets are similar to the linear 
methods for thin sheets except that steps must be taken to account for the 
higher harmonics of the fundamental frequency which are generated, either 
by the material nonlinearity or by other nonlinearity in the experimental 
device such as the acoustic driver in an impedance tube. 

Ingard and Ising (Reference 18) measured the impedance of a single orifice 
at the end of a tube by using a hot wire to measure the velocity in the orifice 
and a pressure pickup in the side of the tube near the orifices. A crude 
Fourier analysis of the signals gave the amplitude of the fundamental and 
the phase of the pressure relative to the velocity. 

Melling (Reference 8 ) has used both the two microphone method and the 
standing wave method to measure nonlinear impedance of perforates. 

Narrow band filters were used to eliminate the higher harmonics such that 
the standard measurements could be made on the fundamental frequency. 

This approach is not reliable, however, because, even though the filters 
remove the higher harmonics from the measurement, these harmonics 
still interact with the fundamental frequency in the material and affect the 
measurement indirectly. If the acoustic driver is given a distorted signal, 
or if the tube is operated near a resonant frequency such that shock Waves 
begin to develop, the measured impedance of the material will show these 
effects even though the signals have been perfectly filtered. Thus, the 
measured Impedance of the material is a function of the entire system 
properties and has no meaning as a material property at all. The solution 
to this problem requires a redefinition of the impedance concept and an 
improved experimental technique. 
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The new impedance concept required for nonlinear impedance was given by 
Zorumski and Parrott in Reference 9 , and is represented by Equation 7. 

This impedance concept used a time domain expression for the material 
property. The acoustic resistance and reactance^ as given by Equation 7, 
depend only on the velocity through the material. An immediate consequence 
of this assumption is that the acoustic resistance must not depend on the 
fundamental frequency or the presence of higher harmonics in an impedance 
tube test. This definition of impedance, if it is correct, is therefore a 
material property Which does not depend on the properties of the measure- 
ment apparatus as does the more conventional definition of nonlinear 
impedance in the frequency domain. Zorumski and Parrott showed that 
their definition does indeed represent a sample of material over a wide 
range of test frequencies and intensity levels, from zero frequency, the 
dc flow resistance, to a 4,000-Hz fundamental frequency test. The 
measurement method used in this test was the two- microphone method with 
the backing cavity adjusted to one-fourth of the fundamental test frequency 
wavelength. The microphone signals were Fourier analyzed and the 
instantaneous pressure differential, velocity through the sample, and 
acceleration through the sample were found to determine the resistance and 
reactance functions in Equation 7. 

The preceding discussion shows that the spectral impedance of a nonlinear 
material depends on the entire frequency spectrum. Wir*^ (Reference 19) 
recognized this intuitively and defined impedance as follows; "The impedance 
at the surface of a nonlinear acoustical material is the ratio of thf pressure to 
the particle velocity, at a particular frequency, and in the presence of a speci- 
fied sound spectrum and a specified turbulence spectrum. " Wirt regarded 
the effect of flow as causing turbulent pressure fluctuations which would 
induce a normal velocity through the material. Although he had no convenient 
way of simulating the flow effect, he did devise a method of measuring imped- 
ance, by the standing wave method, in the presence of an intense broadband 
spectrum. Wirt's impedance tube used three drivers to generate a broadband 
spectrum of specified shape and intensity, while a fourth driver added a pure 
tone which was used for the impedance measurement. One-third octave band 
and 1. 8-Hz bandwidth tracking filters, tuned to the oscillator frequency were 
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used to recover the tone standing ivave pattern. Note that if v (t) in Equation 
25 in an intense broadband spectrum, Equation 25 represents Wirt's defi- 
nition of the impedance. 

2.2.4 Linearization of Impedance 

In this and the previous Section the particle velocity normal to the liner 
surface has appeared in formulae for nonlinear modifications to the liner 
impedance and the same quantity is also used by Rice {Reference 14) to 
relate the flow-induced pressure fluctuations to a flow-induced modification 
to the linear impedance. 


This quantity is not known a priori in practical applications nor is it easily 
measured: the same is true, in fact, for the sound pressure level although 
it is relatively easy to measure and the existence of nonlinear effects makes 
it essential, at least, to estimate the sound pressure level at, or near, the 
liner surface. The relationship between Some sound pressure leVel and the 
particle Velocity is often glossed over in the literature; sometimes the 
plane Wave relation is used (without justification) so that 


* 


V 

cj 


a a 


(26) 


where p* is the pressure amplitude, but of course the correct relation is 
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(by definition). Two problems immediately arise: the pressure p* is the 
complex amplitude actually at the lined surface and is, itself, a function 
of V*. It is not the pressure amplitude in a hard-walled duct and in principle 
it is not known until it can actually be measured with the liner installed. A 
further complication is that the particle velocity, it is argued, at one 
particular frequency should be replaced by a total effective particle velocity 
when using the nonlinear impedance modification formulae, since the 
nonlinear bias is a response to the particle velocity at all frequencies. Thus, 
Rice (Reference 14) and others have defined the total effective particle 
velocity to be 


19 


I 



( 20 ) 


(29) 


(the i subscript denotes a frequency or frequency band and the summation is 
over all the frequency bands which may contribute a significant value for 

In this case p*^^ the pressure amplitude in frequency band i is assumed 
known. Note^ however, that Zi = is calculated using expressions of 
the form 


= [a + b |vg|] - i[g - h |v^|] + i cot(k*d*) (30) 

Then v*^^ and hence v* is calculated using the previous two equations which 
complete^ the first loop; with the new value of the impedance Z^^ can be 
recalculated and the iteration repeated until two successive values of 
are nearly the t.ame. 


In principle, of Course, the pressures p*^ are not known and incorrect trends 
may be obtained from assuming that p*^ is a constant in these iterations. It 
is preferable, following Ingard, (Reference 18) to work in terms of an 
incident sound pressure amplitude, p’!'., and then the particle velocity v’!^. 
is given by 



2p 
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Z. + 1 
1 


(31) 


Twice the incident pressure amplitude, 2p^^, can be loosely identified with 
the pressure amplitude in frequency band i in a given duct system when the 
duct walls are unlined. 
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2. 3 THE EFFECT OF GRAZING FLOW 


2.3.1 Empirical Models 

There is at least one physical effect of sheared mean flow on sound propa- 
gation and attenuation in lined flow ducts which is well established in that 
it has been extensively studied, at least analytically. Physically it is the 
refraction and convection of sound by the mean flow. This effect, when 
included in the propagation equations, yields calculated values for. Say, the 
modal attenuation rates which do not appear to agree with experiment (Kurae 
and Allen, Reference 20), unless a modification is assumed for the liner 
impedance. In addition, from certain experiments, measured liner 
impedances in the presence of grazing flow confirm that a modification 
appears to take place although all experiments to date have been questioned 
in certain respects. While measured attenuation rates do strongly indicate 
some sort of modification this .Tiay be an incorrect interpretation insofar 
as there may be other physical effects which are not accurately represented 
in the propagation equations. However, for the present, this semiempirical 
effect is accepted and expressed as a modification in the form of an 
additional term, to the zero-flow impedance. Incidentally, there is no 
a priori reason why an additive term for the impedance is appropriate. It 
could be equally well that an additional term for the liner admittance is 
appropriate as, for example, in the case where viscous effects on duct 
propagation are represented in the form of an additional admittance 
(Reference 21). 

Kurze and Allen (Reference 20) state that "the reactive component of a lining 
with perforated plate has to be reduced and the resistance had to be chosen 
higher than the liner acoustic resistance of the particular material" to 
obtain good agreement between experimental and theoretical attenuation 
rates. They "assume that the resistance of a nonlinear material increases 
With sound-pressure difference across the lining" but that "the observed 
proportionality factor between the nonlinear influences of flow and sound is 
higher than expected from turbulent pressure fluctuation at a smooth hard 
wall .... The (measured level) is at least 10-15 db lower (however) than 
the equivalent level . . . (but) recent data show that the roughnegs of 
a wall can substantially increase the wall pressure fluctuations. " They 
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conclude that further investigations are required in this subject and that the 
effective resistance should be proportional to the Mach number of the 
mean flow. 


The idea that the turbulent boundary layer may interact with a porous wall, 
in particular a perforate lining, to give rise to the implied effects has been 
developed in at least three ways. 


First, Rice (Reference 14) has argued that the unsteady jets issuing from 

each hole interact with the grazing flow: "pressure oscillations are 

produced within the potential core of the jet due to vortex shedding from the 

grazing flow over this core. " Rice deduces an expression for the meiximum 

value of pressure fluctuations, p , , which is of the form 

g max 


Pg* max 




2 


where 6# is the linear boundary layer thickness. He then shows, with 
certain assumptions, that through the nonlinear effect outlined in the 
previous section the value of the resistance under a typical grazing flow 
velocity, U*, will be 



This linear dependence on the grazing flow mean velocity is in good qualita- 
tive agreement with the direct measurements by, for example, Feder and 
Dean (Reference 22). Rice (Reference 14) agrees with Kurze and Allen 
(Reference 20) that the smooth, hard wall, measured pressure fluctuations 
are too small to explain the grazing flow effect and measurement of pressure 
fluctuations over lined duct surfaces afe urgently required in sufficient 
detail to checkout his proposed model. An identical grazing flow dependence 
is deduced by Eversman (Reference 16) by simply assuming that the grazing 

flow dynamic pressure, p*** is given by 
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where K is a function of the grazing flow boundary layer, face sheet surface 
geometry, and roughness, and is determined empirically. 

The second approach is that developed by Ronneberger (Reference 2^) 
originally for the case where the hole diameter is larger than the boundary 
layer thickness. Based on the concept that the whole boundary layer 
fluctuates in and out of the orifice and on information concerning the 
instability of this motion, Ronneberger deduces that the real part of the 
orifice impedaince is proportional to the grazing flow Velocity, which is 
tojue out very Clearly by his measurements. The dependence on boundary 
layer thickness is not deduced in detail from his theoretical model, but he 
argues that it is of the same form as that deduced by Rice (Reference 14), 
if 6*/dg« 1 andw6*/U*« 1. It should be emphasized that Ronneberger' s 
work is mainly concerned with the case 6*/dg«l and that his postulated 
mechanism differs from that of Rice (Reference 14) in that the latter author 
appeals to the nonlinear response of the perforate whereas Ronneberger 
Views the effect as a flow induced change in the radiation impedance of a 
single orifice. 

Another mechanism has been proposed by Hirata and Itow (Reference 24) 
based on the argument that the kinetic energy of the motion in the vicinity 
of the orifice is "carried away" by the airstream and thus since energy 
"dissipation depends on the diffusion of energy from the resonator into the 
surrounding medium" the dissipation (by the orifice) depends oh the 
(grazing) airflow. They develop simple equations to describe these energy 
related concepts and deduce the same dependence on airflow velocity as the 
previous authors mentioned above. They use their formula for the flow- 
induced resistance in transmission loss calculations which they compare 
with measurements for mean-flow velocities up to 50 m/s and find 
reasonable qualitative agreement. 

Another feature of the flow induced change in acoustic material properties 
is given by Budoff and Zorumski in Reference 25. When there is a tangen- 
tial flow on a side of a thin sheet, it cannot be assumed that the resistance, 
R, for blowing into the stream is the same as for sucking from the stream. 
The flow resistance of perforated sheets was therefore measured with both 
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rxperiments with a single orifice in the presence of grazing flow 
^Reference 22), however, still show a close correlation between acoustic 
resistance and flow resistance. This correlation can be explained from the 
study by Posey and Compton (Reference 26) which shows an equivalent 
resistance can be defined in terms of the resistances for blowing and 
sucking. This equivalent resistance is a weak function of the difference in 
the actual resistances and is roughly 2R^R_/(R^ + R_), where R^ is the 
resistance for sucking and R_ is the resistance for blowing. 

Some unpublished measurements were hnade by Feder (Pratt & Whitney 
Aircraft, East Hartford, CN) of the flow resistance of perforated sheets in 
the presence of tangential flow. Fedef's data were subsequently correlated 
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by Rogers and Hersh (Reference 27) in terms of a discharge coefficient, 
which is related to the specific orifice resistance by 
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For a Sample with <r = 0. 14 1, the resulting resistances become 
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with similar results for other specimens. These results exhibit a dis- 
continuity at the point of flow reversal between blowing and sucking because 
of the singularity in expression for the resistance with blowing R_, but that is 
consistent with the measured results of Budoff and Zorumski (Reference 25). 

It is quite clear from these and other publications that this flow induced 
effect, whatever its physical interpretation, is extremely important in flow 
duct liner design for all practical liner materials, but particularly for 
locally reacting perforate materials, and it is also the least understood. 

2.3.2 Measurement Methods 

A number of methods have been attempted to measure the impedance of an 
acoustic material or an orifice with grazing flow. Before discussing these, 
however, we will briefly mention the far simpler measurement ox flow 
resistance used by Budoff and Zorumski in Reference 25. The experimental 
setup is shown in Figures 6 and 7. The essential part of the method is the 
measurement of the differential in static pressures between the flow duct 
and the side branch and the measurement of steady velocity through the 
sample. Since the pressure differentials may be quite small, it is necessary 
to have a sensor which will indicate differentials down to 10’® atmospheres 
and up to 10 atmospheres. 
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A similar large range of flow rate is required. This was accomplished by- 
using three rotameter flow meters with a suitable pipe and valve arrangement 
as shown on Figure 6. This system also utilizes the same flow meter for 
both blowing and sucking flow rates, so that any asymmetry observed in 
the resistance cannot be due to switching flow meters. 

Feder and Dean (Reference 22) measured the impedance of thin sheets in 
the presence of grazing flow using the standing wave method. The sample 
was mounted in the side of a small flow duct, as in Figure 6, and the side 
branch was used as an impedance tube. The principal difficulties of this 
method were accurate measurement of the standing wave pattern and 
uncertainty about the radiation impedance of the tube, looking into the 
flow duct. 

Ronneberger (Reference 23) used an apparatus similar to Ingard and Ising's 
(Reference 18) to measure the radiation impedance of an orifice looking into 
a shear flow. Some rf his results have been discussed previously in this 
report. We note here, however, that Ronneberger found that the acoustic 
resistance and flow resistance of the orifice were nearly identical and that 
the reactance decreased with increasing tangential flow. 

Armstrong, Beckmeyer, and Olsen (Reference 28) measured the impedance 
of a wave in a flow duct with one treated wall by measviring the axial wave- 
length and decay rate of a duct mode with a two-microphone technique. The 
apparatus is shown in Figure 8. In this method, the shear flow profile is 
measured and the shear flow wave equation is solved to relate the duct liner 
impedance to the measured axial wave number. The advantage of this method 
is that it tests a complete duct liner under conditions simulating its intended 
application. Its disadvantages are that it requires large material samples 
and an elaborate data reduction scheme. (Also, it does not appear to be 
easily extended to the measurement of nonlinear impedance. ) 

Dean (Reference 29) used a one-microphone method (not previously discussed 
here) and the two- microphone method to measure the impedance of duct 
liners in the presence of flow. Dean's results indicated that the two- 
microphone technique was preferable. This method has the advantage of 
being applied to a finished duct liner, like Armstrong's method, and hati the 
added advantage of being able to measure nonlinear impedance. 
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\ K y.i I ‘MMKNniOD MKTllODS FOR MATERIALS 
lUilk ar«‘ nol often used in aiif^raft engine duet applications, 

however, lludr use Is always a possibility so that they should be included 
in any prediction scheme. Little is known about nonlinear behavior of bulk 
materials. The work mentioned here has been concerned with linear 
properties. This linear representation may be acceptable in the aircraft 
engine application because a bulk material in a duct liner weald probably be 
protected by a facing sheet which would cut down on the amplitude of the 
waves in the material. For prediction purposes, therefore, it is 
recommended that: 


Bulk materials are assumed to be linear and described through the equations 


VP + z V =0 

V r-jj u» w 


and 




(35) 

(36) 


where z^ and k^ are experimentally measured complex material parameters 
which are functions of frequency. 

Sheet materials which are essentially rigid, both perforated sheets and 
uniform porous materials, are typically used in aircraft engine as facing 
sheets which are exposed to high intensity sound and to grazing flow. It 
has been demonstrated many times that there is a close correspondence 
between measured acoustic resistance (in the freqvxency domaini and flow 

resistance. Also, time domain measurements have shown this equivalence 
between flow resistance and instantaneous resistance for intense waves 
with fundamental frequencies up to 4,000 Hz. Therefore there is no 
question that the quasi- static assumption is valid over the entire range for 
which measurements can be accurately made. Less is known about the 
reactive part of the material impedance. This term is always small com- 
pared to the resistive part at high amplitudes so that it is extremely difficult 
to measure; however, since it is not a dominant effect there is little motiva- 
tion to measure it and, as an approximation, it may be assumed to be inde- 
dendent of level, that is, the reactive part of the material impedance is linear. 
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It ie known that gra^^ing flow affeccs the resistive part of the impedance of a 
thin sheet. The effects seems to be dlroetly proportional to grassing flow 
velocity (in the free stream) and inversely proportional to boundary layer 
thickness. It also seems that the proportionality constants <sf the flow 
resistance are slightly different for boundary layer blowing than they are 
for boundary layer suction. This difference implies that there is a peculiar 
sort of nonlinearity induced by the grazing flow, even at low sound pressure 
levels. 

The following expression is recommended to represent the acoustic 
properties of all thin sheet materials: 

Ap* = (37) 

where 


PaC 


a '-a 


r[v, M, (^^)J * 


_d_ 
Ca dt* 


(38) 


The nonlinear resistance R is to be measured in a steady flow resistance 

-4 I I -2 

test. The normal velocity v* will have the range 10 < v'^'/c <10 in 


both the blowing (v*/c < 0,M>0) and suction (v*/c >0, M > 0) measurements 

a 


The free stream Mach number shall have the range 0<M<0. 5 and the 
boundary layer displacement thickness shall be chosen to correspond to the 
intended applications of the data. The measured data are to be correlated 
in terms of a discharge coefficient, in terms of which R is given as 
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and the discharge coefficient shall be asymptotic to if |v| «M 

and to a constant when | v | is of the same order as M. In addition, R is 
taken to approach a constant when both M and Jv| are small. The reactive 
part of the impedance, X^, may be measured at a low intensity 
(loo db ^ L s 120 d B) by either the standing wave or the two microphone 


methodi with the material placed in front of a quarter-wave baching cavity. 
The impedance is to be linearized by computing an effective velocity 
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SO that the material impedance is 
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or 


- ix^ (42b) 

If Xg is measured at only one frequency, the expression in Equation 42a can 
be used to calculate at any other frequency by assuming <r, Pe and f* 

are frequency independent (note t = ut /c^). However, it is preferable 
to take measurements of at all frequencies of interest and then to use 
Equation 42b. 


Section 3 
DUCT LINERS 


3. 1 FLAT POINT REACTING LINERS 

Under certain conditions the impedance normal to and at the lined surface 
is a function of only the liner parameters, except for certain local mean 
now properties such as the mean density and the coefficient of viscosity, 
n particular, when it is independent of the properties of the local sound 
field (apart from its frequency) the lined surface is said to be point reacting. 
Physically, the simplest feature of point reaction is that the direction of 

acoustic motion within the liner is independent of the direction of an incident 
plane wave. 


The inepedence normal to and at the surface of a flat, point reacting single- 
layer configuration (backed by an impervious, rigid wall) is of the form. 

= A + i k^/k* cot(k*d*) 

The second term is the impedance of the air-cavity under these conditions 
where d* i, the cavity depth and k* is the wave number of the motion in that 
cavity. In most applications it has been assumed that this motion is in the 
form of a plane wave so that k« = k^ = „/c^. MeUlng (Reference II) has 
i»ted that, due to viscous and other dissipative effects on propagation within 
the honeycomb cavity it is preferable to assign a small Imaginary part to 
this wave number, which could bo obtained from impedance tube or in situ 

measurements, tt only viscous effects are taken into account then (Ref- 
erence 11) k-!* is given by 



(44) 
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wher6 is the effective radius of the honeycomb cell and v is the kinematic 
e 

viscosity of the cavity fluid. The general form of the air-cavity impedance 
expression is for the case when the cavity walls ire normal to the porous 
sheet surface. Appropriate modifications are needed if the cavity has a 
variable cross-sectional area and/or the walls are not normal to the porous 
sheet. Developments in this direction are being pursued by Wirt (Reference 
30) for complicated but point reacting liner configurations which have desir- 
able features such as good low frequency absorption relative to standard 
configurations of the same depth, and wide-band absorption characteristics 
previously thought unattainable with their porous material cavity combination. 


The locally reacting liner may be of the multilayer type, that is, made up 
of a number of porous sheets separated by honeycomb layers of various 
depths. For simplicity we consider the flat single layer type although, in 
principle, the information summarized here can be applied to convex and 
concave, multilayer, locally reacting configurations; general expressions 
for the impedance of those configurations are given by Zorumski (Ref- 
erence 31) in terms of the linear impedance of each material layer and the 
distances between or depths of each layer. Following the analysis of Zwicker 
•»nd Kosten (Reference l7), it can be shown that the impedances of successive 
layers of a duct liner are related by the recurrence formula 
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if the medium within the honeycomb cavities is at ambient conditions. 
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3.2 FLAT EXTENDED REACTION LINERS 

In general, the liner impedance concept is most useful whan the liner is of 
the point reacting type. At the other extreme when the normal impedance 
at the lined surface is, say, a function of the angle of incidence but still 
independent of position, the impedance concept can still be used to some 
advantage, as outlined beloiv. However, in order to incorporate the more 
realistic case when the impedance is a function of position, the impedance 
concept, per se, can be discarded in most extended reaction cases in favor 
of a different approach. That is, Sound propagation through the liner itself 
is included in the description of the wave propagation throughout the duct. 

The boundary condition for the wave prop, ition process is Still in the form 
of an impedance or admittance parameter, but the surface at which it is 
applied is elsewhere such that point reaction can Still be assumed. 

An extended reaction liner is cne in which there may be waves propagating 
parallel to the duct axis, circumferentially around the duct perimeter, or in 
both directions. It is assumed here that the liner is made up of N uniform, 
layers of bulk-type porous material and thin sheets of rigid material such aS 
perforated plates. This liner construction is shown in Figure 9. 
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Whem duct liners are sufficiently thin for curvature effects to be neglected, 
their analysis may be carried out in rectangular coordinates. In this case, 
the wave numbers are related by 

Following the analysis by Zwicker and Kosten (Reference 17), it can be 
shown that the pressures and velocities of successive layers of the duct 
liner are related by the recurrence formula (Equation 45) if is sub- 
Stituted for k|^^ If the liner is partitioned so that waves may not travel 
ajtially, k^^ is set to zero, and if they may not travel circumferentially, 
m is set to zero. Thus, Equation 45 provides a means of calculating the 
surface impedance of uniform, multi-layered extended reaction liners. 

Through Equation 46 the surface impedance of such a duct liner may be a 
function of the circumferential and/or axial wave numbers. This generally 
presents no difficulty as far as the circumferential wave number is concerned, 
since the acoustic equations in aixisymmetric ducts may be solved separately 
for each circumferential wave. However, the impedance Equation 45 must 
be evaluated simultaneously with the appropriate duct propagation equations 
and k^^ will be Variable during the solution process. 

3. 3 recommended METHODS FOR LINERS 

An iteration process must be uSed to compute liner properties since the 
impedance depends on the acoustic spectrum through Equation 40. Here we 
assume that a one-third octave band spectrum of p^ is known at the linear 
surface. The first step is to compute the resistances of the material layers 
from some estimate of the corresponding velocities. This information is 
then used to re-estimate the velocities. The process may be started by 
assuming all velocities are small so that the material impedances or resist- 
ances are independent of velocity. Thus, for a liner with N layers, we start 
with the velocity vector for each frequency or frequency band i 
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Thdn Equation 40 id 



used to computa an effective velocity for each layer 
n 1/21 



(48) 


and Equationd 41 and 42 are used to compute an impedance vector, at 

each frequency, ivhere ’ * 
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( ^ ) 


1 e e j 


(49) 


The condition p. = p^^ and vj^^ = 0 are used with sets of Equations 45 
to give 


D 0] 




Pi 

[T?J.O 
• • • 

◄ 

(N) 
v; ' 
1 

• 

• 

^ 

0 

* 

> 

• 

M [Ti] 

[o 9 


p!°> 

vP> 

^ » J 

1 

0 

^ J 


where the identity matrices in Equation 50 are order (2 x 2) and the square 
matrices (J„_ 2 ) are the products shown in Equation 45. The solution to 
Equation 50 at each frequency gives a new set of velocity vectors which may 
be used to repeat the iteration starting with Equation 47. The iteration 
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stops when the effective velocity vector ceases to change by a significant 
amount with successive Iterations (or, alternatively, the re sistance vector 
as described in Section 4). 


3.4 EXAMPLE LINER 

Reference 32 describes the inlet liner for a typical ground lest nacelle. 

The liner is a perforated plate, 0. 5-mm thick, with 1. 27-mm diameter 
holes spaced to achieve an open area <r = 0. 068. The backing is honeycomb 
with a 12.2-mm depth. The needed flow resistance and effective density 
function for the plate is not known, however they may be estimated. For 
flow into the liner the orifice discharge coefficient is ass umed to be 
0. 8 /vTm (Reference 27) and for outflow it is taken as >/|v|/M. The corre- 
sponding orifice resistances are M/1. 28 and M/2. 0 and the acoustic resist- 
ances (based on average velocity through the plate) are = M/0. 087 for 
inflow and R ^ M/0. 136 for outflow. The effective density is taken to be the 
air density so that = 1. The inlet mach number is 0. 58 so that the effec- 
tive resistance is R^ = 5. 20. The reactance in Equation 42 is negligible so 
that Equation 43 gives an impedance of 

Z = 5. 20 + i cot (2. 253X10"'^f), 
which is 5. 20 + i 4. 36 at f = 1 kHz. 
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Section 4 

RECOMMENDED PREDICTION METHODS 


4. 1 SINGLE-LAYER LINER CONFIGURATIONS 

The calculation method falls into two parts. The first is concerned with the 
cavity impedance, which is assumed to be independent of sound pressure 
level and grazing flow effects. The second involves an iteration for the sheet 
material acoustic resistance which is dependent on those effects and which 
also utilizes the calculated cavity impedance. Input data for both calculations 
is summarized first, under the data base heading. 


It should be noted that symbols like Z^, p.^ and v. have been replaced by 
Z(f.), Pq(I^) and that is, v^, for example, at frequency band i is now 

particle velocity v(L) at the frequency band j for which the center frequency 
is fj. This not only serves to provide an explicit reference to f. (and pre- 
vents confusion between i = \ -1 and the frequency number), but it dis- 
tinguishes the parameters used here from their counterparts elsewhere in 
the report which are often single frequency or spectral (Fourier) quantities. 
Also is should be noted that, in general, the cotangent function in what follows 
and the parameters k»:=(L), z*(L), Z^(f.), k^^; AZ(L) and v(f.) are 

complex quantities. 


Data Base 


1 . 

2 , 

3. 

4. 


Nf. 


Hz 


Center-band frequencies L, j = 1, 2 
Speed of sound in duct medium, c_, m/sec 

cl 

Characteristic impedance of duct medium, p c , 

21 3 . 


MKS RAYLS 


Sound pressure level spectrum in hard-walled duct, SPLj^j (L), 


dbre: p^.^^ 


= 2 X 10“® N/m^ 


5. Incident (normalized) sound pressure spectrum: 




xPref' P 


a a 


with SPL,j(i.) from (4), p from (2), (3) 
1 j j a a 


39 


6. 

7 . 

8 . 


9 . 


10 . 


11 . 


R versus v material D. C, flow characteristic (a1 tin* n quirofl 
tangential Mach number) 

Material specific acoustic reactance at each frequfiu'v, X (f.) 

' t * 1 

Wave number and distributed impedance of cavity filler (if any) at 
each frequency, k-:«(r), rad/m and z'^'ifj), MKS RAYLS/m 
Cavity depth, d'i*, m if simple cavity configuration (see (12) below). 
Honeycomb effective radius, r *, m if simple, small- radius 
honeycomb forms air-filled^ cavity. 

Kinematic viscosity of cavity medium, v. vn /b&c if simple, small- 
radius honeycomb forms air-filled^ cavity. 


The specific acoustic cavity impedance, Z (f.), is determined at 


Calculation of Cavity Impedance 

12 , 

each frequency with one of the methods outlined below, as 
appropriate. 

a. Simple air-filled cavity; 


Z^(f.) := i|k^/k*(f.)j cot jk*(f.)d=:j 
where 

k*(f.) = k = 2rf./c^ 

J a J a 


or for a "small radius" honeycomb cavity 


27rf. 

k*(f.) ^ + 

J "a 


(-fjV) 


1/2 




b. Simple cavity, filled with bulk material : 

Z^j(fj) = -2*(r)/|p^c^ |k*(fj)d>:j 

c. Complex, point reacting cavity geometries, curve d liner 
configurations, etc. — use measured values for Z (f .) or 
values calculated according to referenced theoretical models. 

^i. e. , the cavity medium is acoustically identical to the duct medium. 
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d. Extended reaction (nonpartitioned) liner configuratione — the 
impedance calculation method must be carried out in conjunc* 
tion with the solution of the duct propagation equation, which 
provides values for m/b. In this case Z^(fj) is 

calculated with 

(f.l)]cot {kj„ «j) 


where 


''Jm <9 /’'a' 


= {(k* (tj)/kj2 . (f.) . 6^(m/b)^ 


1/2 


(k^=2rf./cj 

and 6.6 are unity in the absence of partitioning, but 6=0 
z rti z 

if the liner is axially partitioned or 6^ = 0 if circumferentially 
partitioned. 

Liner Impedance Calculation: Iteration for Sheet Material Acoustic 
Resistance 

13. Choose an initial value for the sheet material specific acoustic 
resistance, e. g . , 


R„ = 1.0 

e 

14. Form a value for the complex material specific acoustic impedance 
at each frequency 

AZ(f.) = R,, - i X^(f.) 

J e e j 

with R from (13) or (18) and X (f.) from (7). 
e e j 
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^5. Foytn a value for the complex liner specific acoustic impedance at 
each frequency 

Z(f.) = AZ(f^) • 7,^(f.) 

with AZ(fj) from (14) and Z^(fy from (12). 

16. Calculate the complex effective particle velocity Mach number at 
each frequency 

V (fj) = 2p^ 

Pq (^i) from (5) and Z(f.) from (15). 

J 

1 1 . Calculate the total effective particle velocity Mach number 



j = 1 


with V (fj) from (16). 

13. Read off from R[v], the material D. C. How characteristic 
(from (6)), the values R^, R where 



and calculate the effective specific acoustic tesistance 
Re = 2R^R_/(R^ + R^). 

If :his value of R^ differs from its previous value (or if it differs from that 
ass’-cned in (13) by less than a specified tclc’-mce then the iteration can be 
terminated and the required specific acoustic liner impedance is the last 


value calculated at (15), Otherwise the procedure must be repeated starting 
at (14) with the new value of or one based upon this and previous values 
(for example by linear extrapolation). 

4. 2 MULTILAYER LINER CONFIGURATIONS 

The calculation method falls into two parts. The first is concerned with the 
cavity impedance matrices^ which are assumed to be independent of sound 
pressure level and grazing flow effects. The second involves an iteration 
for the sheet material resistances which are dependent on the former effec*, 
while the surface layer alone is influenced by grazing flow effects. Input 
data for both calculations is summarized first, under the data base heading. 

The reader is referred to the comments on notation, etc. , in the introcuctirr. 
to Section 4. 1. 


Data Base 

1. Center-band frequencies f., j = 1, N^, Hz. 

2. Speed of sound in duct medium, c , m/sec. 

cl 

3. Characteristic impedance of duct medium, p c , MKS RAYLS. 

3l SL 

4. Sound pressure level spectrum in hard-walled duct SPL^^ (f.), 
db re; Pj.g£ = 2 x 10"^ N/m^. 

o. Incident (normalized) sound pressure spectrum: 





10 


SPL^(f .) /20| 




2 

a 


with SPL„ (f.) from (4), p c^ from (2), (3). 

n J ci SL 

6. a. The number of material layers, N 

b. R' ' versus v material D. C. flow characteristic for each 

material layer (at the required tangential Mach number for the 
surface layer, n = N) 

7. Material specific acoustic reactance at each frequency for each 

layer, (f.) 

8. Distributed impedance and wavenumber of cavity filler 'if any) at 

each frequency for each layer, (f.), rad/m and z**”' (f.), 

MKS-RAYLS/m. 
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f 


9. Cavity depth of each layer, m, If simple cavity configurations 

(see (12) below). 

*«» / V 

10, Honeycomb effective radius, r , m, for each layer, if simple, 

® t 

small-radius honeycomb forms air-filled cavities. 

11. Kinematic viscosity of cavity medium, v, m^/sec, if simple, 
small-radius honeycomb forms air-filled^ cavities. 


Calculation of Cavity Impedance Matrices 

12. The cavity impedance matrices are defined by 


a(h) 


(fi) 


(n) 


cos k 


sin k 


(n) 



sin k 
cos k 


in)' 

(n) 


where 

^(n) ^ ^=:=(n)^ n = 1, N 

^(n) ^ k''^*'^f.)j . n = 1, N 


The wave numbers k and distributed impedances z for each 

J » 

layer, at each frequency are calculated with one or more of the methods 

outlined below, as appropriate. 

a. Simple air-filled cavity: 


k*^”^fj) = 2-f./c^ 


or for a small-radius honeycomb cavity 

, 1/2 


k '-"(fj) — 


2rrf. (rf.v) 

+ i 


/f \ 1 X ; I 


^Hti) 
e a 


and 

= -i(2rf./cj p^c^ 


^i. e. , if the medium in all the cavities is identical acoustically to the duct 
medium. 


.1 • , ; ^♦*...5 

•V H -V, 
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I 


t. / \ p'* / 

b. gimple cavity filled with bulk mater ialj k‘ and z' '(£j) 

£ollow directly £rom (8). 

c. Complex, locally reacting cavity geometrioe, curved liner a 

con£ieuration8. etc. the matricoB muat be derived 

directly from measurement and/or theoretical models. 

d. Extende d reaction (nonpartitioned) liner configurations — the 

impedance calculation method must be carried out in conjunc- 
tion with the solution of the duct propagation equation, which 
provides values for k and m/b. In this case a can be 

obtained from a or b above while k ' (£^) is denoted by 

^rm^ (fj) and is given by (for all layers) 


Mn) 


.(n) 


k '“'(f;)/k^ = k';;'(f.) 
rm j a rm j 


(k"^‘"^f.)/k Y - kf^(f.) - 5^(m/b)^ 

y ' j a/ z zm j m 


1/2 


n = 1, N 


and k’ ^”^(fj) can be determined from a or b above; 6^, 6^ are 
unity in the absence of partitioning, but = 0 if the liner is 
axially partitioned or 6^^ = 0 if circumferentially partitioned. 


Liner Impe dance Calculation: Iteration for Sheet Material Acoustic 
Resistances 

13. Choose an initial value for the material specific acoustic resistance 


for each sheet e. g. , 





1 \ 


r(N) 

e 


1.0 

K"1 - ■ 

j,(N.l) 

e 

► r i 

1.0 

• 


r('d 

V^e ^ 


• 

^l.O; 
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14. Form a value for the complex material specific acoustic impedance 
for each sheet at each frequency 

- i X^"^f.)) 

with from (13) or (18) and jx^g"Nf.) | from (7), 

15. Calculate the complex liner specific acoustic impedance at each 
frequency by repeated application of the following equation, start- 
ing with (f.) = 0, p^®^ (f.) = 1.0 

J J 


J 





1 AZ^®V.) 

J 


A^"\f.) 

J 


1 


with (f.) from (12). (The resulting values (f ), v^^^ (f.) 

J (o ) ' ’ 

correspond to an arbitrary value of p ’ (f.) and are denoted by tilde, ) 

The liner specific acoustic impedance is given by 


Z(f.) = p^^'*(fy/v'^^f.) 


In. a. Calculate the actual, complex particle velocity Mach number at 
each frequency for the surface layer (n = N) 

(f.) = 2p^(fy/{z(f.) + ij 

with p (f.) from (5) and Z(f.) from (15). 

^ J J 

b. Calculate the actual, complex particle velocity Mach number 
for all the remaining layers at each frequencv 



with v^^\f.) from (15) and v^^\f.) from (16a). 
J J 
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17. Calculate the total effective particle velocity Math number for each 
layer 




1 

v<“',ry 

2 

1/2 





• 



with Iv'”^ (f.)l from (16b). 

J ) /j^\ 

18. Read off from R' ’ [v] (the material D. C. flow characteristics 

from (6), the sheet material specific acoustic resistances corres- 
ponding to positive and negative values of the effective velocity 
Mach numbers from (l7): 


j,(n) 

II 

ja 

1 1 

+ 

< 

0£ 

j^(n) 

= l4-v<">i! 


1 L ® J 


and calculate the effective specific acoustic resistances with 

If each value R^ , n = 1, N differs from its previous value (or if 
each differs from that assigned in (13) by less than a specified 
tolerance then the iteration can be terminated and the required 
liner specific acoustic impedance is the last value calculated at (15). 
Otherwise the procedure must be repeated starting at (14) with the 
new values or a vector based on this and preceeding vectors 

(for example, by linear extrapolation). 
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Section 5 

CONCLUDING REMARKS 


A natural result of a review is the identification of unanswered questions. 

Many such questions could be raised as to why a material has a certain 
acoustic property, but in this paper this area has been neglected, assuming 
that the material property would be measured and then used to compute the 
liner property. This is a valid approach and is recommended as the method 
to be used in the future; however, it is recognized that this approach requires 
both the development and standardization of precise measurement techniques. 

In view of the importance of flow on material properties, it is recommended 
that precision methods be developed for the measurement of the flow resistance 
of materials in the presence of flow, including the boundary layer parameters 
as variables. It is also recommended that methods of measuring the nonlinear 
acoustic impedance of materials in presence of flow be developed. This impe- 
dance measurement must be done using time-domain techniques to see if the 
correspondence between flow resistance and acoustic resistance is retained 
valid in the presence of tangential flow. 



Very little is known about the way a broadband spectrum affects the apparent 
(linearized) properties of a material. In the present paper, an rms effective 
velocity was recommended because it is realized that complex spectra, or 
real-time information is not generally available. It is Worthwhile exploring 
some of these questions by looking at the material and liner properties the 
time domain, where the nonlinear effects are easier to characterize but the 
geometric effects, such as liner depth, are more complicated. 


The extcnsioft of the properties of linear locally reacting liners to the case of 
extended reaction is deceptively simple. It should be emphasized that this is 
not possible in the case of nonlinear liners. Further research is needed on 
nonlinear extended reaction liners; however, it is beyond the scope of the 
present paper to discuss this in detail. 
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